Theoretical analysis for the apparent discrepancy between pp and pp data in charged 
particle forward-backward multiplicity correlations 
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The strength of charged particle forward-backward multiphcity correlation in p + p and p + p 
coUisions at ^5=200 GeV is studied by PYTHIA 6.4 and compared to the UA5 and STAR data 
correspondingly. It is turned out that a factor of 3-4 apparent discrepancy between UA5 and STAR 
data can be attributed to the differences in detector acceptances and observing bin interval in both 
experiments. A mixed event method is introduced and used to calculate the statistical correlation 
strength and the dynamical correlation strengths stemming from the charge conservation, four- 
momentum conservation, and decay, respectively. It seems that the statistical correlation is much 
larger than dynamical one and the charge conservation, four-momentum conservation and decay 
may account for most part of the dynamical correlation. In addition, we have also calculated the 
correlation strength by fitting the charged particle multiplicity distribution from PYTHIA to the 
Negative Binomial Distribution and found that the result agrees well with the correlation strength 
calculated by mixed events. 
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Fluctuations and correlations are important observ- 
ables investigating the properties of thermodynamic sys- 
tem and are critical tools revealing the mechanism of 
particle production and the formation of quark-gluon- 
plasma in relativistic heavy ion collisions Several 
thermodynamic quantities and the produced particle dis- 
tributions show varying fluctuation patterns when sys- 
tem undergoes phase transition. Such as the large energy 
density fluctuation is expected in the first order phase 
transition and a second order phase transition may relate 
to a divergence in specific heat. The event-by-event fluc- 
tuation pattern in average transverse momentum may 
significantly change around a critical point, etc. 

The experimental study of correlation and fluctuation 
becomes a hot topic in relativistic heavy ion collisions 
with the availability of high multiplicity event-by-event 
measurements at the CERN-SPS and BNL-RHIC exper- 
iments. There have accumulated an abundant experi- 
ment data d, i, i, i, 0, i, H mini where arisen new 
physics are urgent to be studied. 

Recently STAR collaboration has measured the 
strength of charged particle forward-backward multiplic- 
ity correlation, h (deflned later), yh p + p collision at 
^/s=2m GeV [il [H. It is 3-4 times smaller than the 
one measured by UA5 vci p -\- p collision at the same en- 
ergy apparently liBl- In this paper, the PYTHIA 
6.4 [i^ is employed to analyze both STAR p + p and 
the UA5 p + p data. It is turned out that the above 
apparent discrepancy is because of the differences in de- 
tector acceptances and the interval of pseudo-rapidity 



bin in both experiments. In addition, a mixed event 
method is proposed and used to calculate the statistical 
correlation and dynamical correlations stemming from 
the charge and four-momentum conservations and the 
decay of unstable particles individually. We also flt the 
particle multiplicity distribution from PYTHIA to the 
Negative Binomial Distribution (NBD) and calculate the 
strength of charged particle forward-backward multiplic- 
ity correlation which agrees well with the one calculated 
by mixed events. 

Following Refs [H, the strength of charged particle 
forward-backward multiplicity correlation, 6, is defined 



h = 



(nfTit,) - {nf)(ni,) _ cov{nf,ni,) 



(("/ - W»K - ("b>)> 



(1) 
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where n/ and rib are, respectively, the number of charged 
particles in forward and backward pseudo-rapidity bins 
(A77) defined relatively and symmetrically to a given 
pseudo-rapidity rj. The (n/), for instance, refers to the 
mean value of n/ and the cov(nf,nh) and var(nf) are, 
respectively, the forward-backward multiplicity covari- 
ance and forward multiplicity variance. If there is no 
correlation between forward and backward multiplicity, 
then (nfTib) = {nf){ni,) and 6 = 0. Thus 6 is a measure 
of the strength of forward-backward multiplicity correla- 
tion. As the denominator in Eq. [T]is positive, if both Uf 
and Uf, are, respectively, larger or smaller than (rif) and 
(rib) simultaneously the correlation is positive, negative 
otherwise. 
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Although both STAR and UA5 experiments ^ \m, 
[T^ . [isj measure the charged particle multiplicity in Non- 
Single-Diffractive (NSD) p + p and p + p collisions at 
•y/i=200 GeV, the detector acceptances are different 
from each other. In UA5 ex per iment they are pt > 
GeV/c and 0.0 < |r;| < 4 [lljll^.ut pr > 0.15 GeV/c 
and 0.0 < |r/| < 1.0 in STAR[iljll|. Meanwhile, the ob- 
served interval of forward- backward pseudo-rapidity bin 
is Ar/=1.0 in UA5 experiment rather than 0.2 in STAR. 
We shall show that those differences are the origin of 
apparent discrepancy between STAR and UA5 data in 
the strength of charged particle forward-backward mul- 
tiplicity correlation. 

The UA5 data of energy dependence of the correla- 
tion strength have been studied by Dual Parton Model 
[13 and the statistical model |ll| and the STAR data 
of forward-backward charged particle multiplicity covari- 
ance have been studied recently in fl^. However, the ap- 
parent discrepancy in correlation strength between UA5 
and STAR data is not investigated yet. In this paper 
PYTHIA 6.4 fie] is employed to study that. Since we 
are not aim to reproduce the experimental data but to 
study the physics, we do not adjust the model parame- 
ters and default values are used in all the calculations. 

The comparison of experimental strength of the 
charged particle forward-backward multiplicity correla- 
tion to the corresponding theoretical results is given in 
Figure [T] where the upper panel is for UA5 p+p collision 
and the lower panel for STAR p + p collision at y/s=200 
GeV. One sees in this figure that the theoretical results 
are not so far apart from the experimental data for both 
the p + p and p + p collisions. The theoretical correlation 
strength in p + p collision is also a factor of 3-4 larger 
than the one in p + p collision, especially. 

In the upper panel of Figure [2] the full squares are 
the PYTHIA results for p + p collision with same detec- 
tor acceptances and rj bin interval as in UA5 experiment, 
whereas the open triangles are the PYTHIA results with 
PT > 0.15 instead of pT > GeV/c. The open trian- 
gles are monotonously below the full squares. Middle 
panel of Figure [2] shows the PYTHIA results calculated 
at the same detector acceptances as UA5 but with var- 
ied pseudo-rapidity bin intervals: A77 =1.0 (full squares), 
0.5 (open circles), and 0.2 (open triangles-down), respec- 
tively. Here one knows that the correlation strength, 
b, declines dramatically with the decreasing of pseudo- 
rapidity bin interval. The PYTHIA results calculated for 
both the p + p (full squares) and p + p (open triangles- 
up) collisions at STAR detector acceptances and pseudo- 
rapidity bin interval are given in lower panel of Figure 
[21 We see in this panel that a factor of 3-4 apparent 
discrepancy nearly disappears if both STAR and UA5 
experiments are performed at the same detector accep- 
tances and pseudo-rapidity bin interval. 

It is very hard to separate the statistical and dynam- 
ical correlations (fluctuations) from the measured corre- 
lations (fluctuations) [2]. We introduce a mixed event 
method based on the real (PYTHIA) events (150 thou- 
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FIG. 1: The strength b of the charged particle forward- 
backward multiplicity correlation in p + p (upper panel) and 
p + p (lower panel) collisions at ^5=200 GeV. The experi- 
mental data are taken from [l^ and [l^ . respectively. 



sand events, for instance). The mixed events are gen- 
erated one by one according to the real events. We as- 
sume flrst that the particle multiplicity V in a mixed 
event is the same as the corresponding one in real events. 
However, the N particles in a mixed event are sampled 
randomly from the particle reservoir formed by all par- 
ticles in the real events. As the particles in a mixed 
event are separately and randomly taken from different 
real events, there is not any dynamical relevance among 
them. Thus the correlation strength, b, calculated by 
the mixed events is reasonably to be identified as the 
statistical correlation. Of course, we can also gener- 
ate the mixed event with individual constraint, such as 
charge conservation, four-momentum conservation, and 
decay, etc. The corresponding correlation will be indi- 
cated by "statistical plus charge dynamical correlations" , 
the "statistical plus four-momentum dynamical correla- 
tions", and the "statistical plus decay correlations", etc., 
respectively. 

The strength of charged particle forward-backward 
multiplicity correlation is calculated individually from 
the real events, mixed events, mixed events with charge 
conservation, and the mixed events with charge and four- 
momentum conservations. They are given in the upper 
panel of Figure [3] by the full squares (indicated as total 
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FIG. 2: Correlation strength 6: upper panel for p+p collision 
calculated at different pr cut; middle panel calculated at var- 
ied intervals of rj bin; lower panel for p + p and p + p collisions 
calculated at STAR detector acceptances and pseudo-rapidity 
bin interval. The reaction energies are all y^=200 GeV. 



correlation), open circles (statistical correlation), open 
triangles-up (statistical plus charge dynamical correla- 
tions) , and the open triangles-down (statistical plus the 
charge and four- momentum dynamical correlations), re- 
spectively. In lower panel of Figure [3] are given the total 
dynamical correlation strength (full squares), the charge 
dynamical one (open triangles-up), the four- momentum 
dynamical one (open triangles-down), and the decay dy- 
namical correlation (open circles). The former three are 
extracted from upper panel by subtracting, respectively, 
the statistical correlation from total correlation, the sta- 



tistical one from "statistical plus charge dynamical" one, 
the "statistical plus charge dynamical" one from "statis- 
tical plus charge and four-momentum dynamical" one. It 
has to mention that in the real events generated above 
the decay of unstable hadrons is allowed. In order to 
calculate the decay dynamical correlation (open circles 
in lower panel of Figure [3]), we have first to generate 
the real events without decay of unstable hadrons and 
the corresponding mixed events. One subtracts the b 
calculated by mixed events with (without) decay from 
the b calculated by real events with (without) decay, 
one has the dynamical correlation with (without) decay 
then. The decay dynamical correlation in lower panel 
of Figure [3] is just resulted by subtracting the "dynami- 
cal correlation with decay" from the one without decay. 
One knows here that the dynamical correlations stem- 
ming from charge, four-momentum, and the decay may 
account for the most part of the total dynamical corre- 
lation. 
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FIG. 3: Statistical and dynamical correlation strengths. See 
text for detail. 



Many experiments have indicated that the parti- 
cle multiplicity distribution in hadron-hadron and the 
nucleus-nucleus collisions is well described by the Nega- 
tive Binomial Distribution (NBD) [13, [Ml [gl, [H . For 
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an integer n the NBD reads 



n + k -1 
k - 1 



(1 + Ai/fc)"+'''' 



(2) 



where fi = (n) is a parameter, P(n; /z, fc) is normahzed in 
< n < oo, and k is another parameter responsible for 
the shape of the distribution. If fc is a real the binomial 
coefHcient in Eq. [2]is fc(fc + 1) • • • (fc + n - 1) /n\. In NBD 
the variance (ct^) and mean (/i) is related to k by 



(3) 



An important property of NBD is that if particle multi- 
plicity n is NBD in whole phase space and the particle 
has unified probability, p, in a partial phase space (such 
as in a ?7 bin here) then particle multiplicity distribu- 
tion in this partial phase space is also NBD with same 
parameter k and the mean is equal to fj,p [23 |. That is 
obviously based on the assumption that the particles are 
independent with each other, i.e. there is no dynamical 
correlation among them. 

Since forward and backward pseudo-rapidity bins are 
symmetry relative to the investigated pseudo-rapidity 
[r]) and have same width so (nf) = (rib) and var(nf) = 
var{nb) in NBD. Using the statistic formula [1^ 

var{X ± r) = var{X) + var{Y) ± 2cov{X, Y) (4) 



the b can be written as 



6 = 



cov{nf,nb) var{nf + rifi) ^ 2var{nf) 



varijif) 2var{nf) 
Substitute Eq. [3] into Eq. [5] one has 



(n/) + k 



(5) 



(6) 



We know that the NBD becomes a Poisson distribution 
in the limit k ^ oo, so the correlation strength, 6, is zero 
in Poisson distribution. If the charged particle multiplic- 
ity distribution in real events with decay assumption is 
fitted by NBD, the parameter k is obtained. As the real 
events are generated in NSD (Non-Single-Diffractive) in- 
deed, the charged particle multiplicity distribution is not 
perfect NBD, therefore the above fit is not so sensitive to 
the k values within 6-7. If NBD with k—Q.& is assumed 
for the charged particle multiplicity distribution, the cor- 
responding h can be calculated by Eq. [6]because (n/) can 
be approximated by dNch/dt] in real event. Those b are 
shown in Fig. |4]by open triangles. In this figure the full 
squares are calculated by the mixed events with decay 
assumption (i.e. the open circles in upper panel of Fig. 
121) and the open circles are the charged particle pseudo- 
rapidity distribution in real events with decay (in draw- 
ing dNch/drj the abscissa is identified as ij and scaled 
by 2). The results of NBD agree well with the results 
calculated by the mixed events with decay assumption. 



it proves again that we are reasonable identifying the b 
calculated by mixed events as the statistical correlation 
strength. Comparing the full squares and open triangles 
to the open circles one knows that the statistical corre- 
lation strength may have shape similar to the charged 
particle pseudo-rapidity distribution. 

Recently, the forward-backward multiplicity covari- 
ance in p+p collision at ^/s=200 GeV has been studied in 
Ref . [101 . They assumed the back-to-back partonic scat- 
tering is the origin of hadronic correlation, related that 
partonic scattering angles to a Gaussian like hadroniza- 
tion function, and derived the forward-backward multi- 
plicity covariance. Without more dynamical input their 
results are well comparing with STAR data [l^l . There- 
fore they conclude that the correlation length might have 
no fundamental significance. We plan to investigate the 
partonic origin of forward-backward multiplicity corre- 
lation by transport model in next study. 
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FIG. 4: Charged particle forward-backward multiplicity 
statistical correlation strength b and the charged particle 
pseudo-rapidity distribution in p+p collision at ^5=200 GeV. 

In summary, we have calculated the strength of 
charged particle forward-backward multiplicity correla- 
tion in p + p and p+p collisions at -^5=200 GeV by 
PYTHIA 6.4J1II and compared with UA5 data [3 and 
STAR data [I4I, respectively. It is turned out that a 
factor of 3-4 apparent discrepancy between UA5 and 
STAR data can be attributed to the differences in de- 
tector acceptances and the interval of observed i] bin 
in both experiments. A mixed event method is intro- 
duced and used to calculate the statistical correlation 
strength and the individual dynamical correlations stem- 
ming from charge conservation, four-momentum conser- 
vation, and the decay. It seems that the statistical cor- 
relation is much larger than the dynamical one, and the 
charge, four-momentum, and decay may account for the 
main part of the dynamical correlation. The NBD b re- 
sults agree well with the ones calculated by mixed events 
proves again that one is reasonable to identify the cor- 
relation in mixed events as a statistical correlation. 
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